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INTELLIGENT ROBUST CONTROL SYSTEM FOR MOTORCYCLE USING SOFT 

COMPUTING OPTIMIZER 

5 

Background 

Field of the Invention 

The present invention relates to electronically-controlled motorcycle steering and 
maneuverability systems based on soft computing. 

10 Description of the Related Art 

Feedback control systems are widely used to maintain the output of a dynamic system at 
a desired value in spite of external disturbances that would displace it from the desired value. 
For example, a household space-heating furnace, controlled by a thermostat, is an example of 
a feedback control system. The thermostat continuously measures the air temperature inside 

15 the house, and when the temperature falls below a desired minimum temperature the 
thermostat turns the furnace on. When the interior temperature reaches the desired minimum 
temperature, the thermostat turns the furnace off. The thermostat-furnace system maintains the 
household temperature at a substantially constant value in spite of external disturbances such 
as a drop in the outside temperature. Similar types of feedback controls are used in many 

20 applications. 

A central component in a feedback control system is a controlled object, a machine or a 
process that can be defined as a "plant", whose output variable is to be controlled. In the 
above example, the "plant" is the house, the output variable is the interior air temperature in 
the house and the disturbance is the flow of heat (dispersion) through the walls of the house. 

25 The plant is controlled by a control system. In the above example, the control system is the 
thermostat in combination with the furnace. The thermostat-furnace system uses simple on-off 
feedback control system to maintain the temperature of the house. In many control 
environments, such as motor shaft position or motor speed control systems, simple on-off 
feedback control is insufficient. More advanced control systems rely on combinations of 

30 proportional feedback control, integral feedback control, and derivative feedback control. A 
feedback control based on a sum of proportional, plus integral, plus derivative feedback is 



often referred as a linear control. Similarly, Proportional feedback control is often referred to 
as P control, and Proportional plus Derivative feedback is often referred to as P(D) control, 
and Proportional plus Integral feedback is referred as P(I) control. 

A linear control system (e.g., P, P(I), P(D), P(I)D, etc.) is based on a linear model of 
5 the plant. In classical control systems, a linear model is obtained in the form of ordinary 
differential equations. The plant is assumed to be relatively linear, time invariant, and stable. 
However, many real-world plants are time varying, highly non-linear, and unstable. For 
example, the dynamic model may contain parameters (e.g., masses, inductance, aerodynamics 
coefficients, etc.), which are either only approximately known or depend on a changing 

10 environment. If the parameter variation is small and the dynamic model is stable, then the 
linear controller may be satisfactory. However, if the parameter variation is large or if the 
dynamic model is unstable, then it is common to add Adaptive or Intelligent (AI) control 
functions to the linear control system. 

AI control systems use an optimizer, typically a non-linear optimizer, to program the 

15 operation of the linear controller and thereby improve the overall operation of the control 
system. 

Classical advanced control theory is based on the assumption that all controlled "plants" 
can be approximated as linear systems near equilibrium points. Unfortunately, this assumption 
is rarely true in the real world. Most plants are highly nonlinear, and often do not have simple 
20 control algorithms. In order to meet these needs for a nonlinear control, systems have been 
developed that use Soft Computing (SC) concepts such Fuzzy Neural Networks (FNN), Fuzzy 
Controllers (FC), and the like. By these techniques, the control system evolves (changes) in 
time to adapt itself to changes that may occur in the controlled "plant" and/or in the operating 
environment. 

25 Control systems based on SC typically use a Knowledge Base (KB) to contain the 

knowledge of the FC system. The KB typically has many rules that describe how the SC 
determines control parameters during operation. Thus, the performance of an SC controller 
depends on the quality of the KB and the knowledge represented by the KB. Increasing the 
number of rules in the KB generally increases (very often with redundancy) the knowledge 

30 represented by the KB but at a cost of more storage and more computational complexity. 
Thus, design of a SC system typically involves tradeoffs regarding the size of the KB, the 
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number of rules, the types of rules, etc. Unfortunately, the prior art methods for selecting KB 
parameters such as the number and types of rules are based on ad hoc procedures using 
intuition and trial-and-error approaches. 

Steering and/or maneuverability control of a motorcycle using soft computing is 
5 particularly difficult because of the difficulty in obtaining a sufficiently optimal knowledge 
base (KB). If the KB does not coontain enough knowledge about the dynamics of the 
motorcycle-rider system, then the soft computing controller will not be able to steer and/or 
maneuver the motorcycle in a satisfactory manner. At one level, increasing the knowledge 
conatined in the KB generally produces an increase in the size of the KB. A large KB is 
10 difficult to store and requires relatively large amounts of computational reasources in the soft 
computing controller. What is missing from the prior art is a system and method for 
controlling a motorcycle using a reduced-size KB that provides sufficient knowledge to 
provide good control. 

Summary 

15 The present invention solves these and other problems by providing a SC optimizer for 

designing a globally-optimized KB to be used in a SC system for an electronically-controlled 
motorcycle. In one embodiment, the SC optimizer includes a fuzzy inference engine. In one 
embodiment, the fuzzy inference engine includes a Fuzzy Neural Network (FNN). In one 
embodiment, the SC Optimizer provides Fuzzy Inference System (FIS) structure selection, 

20 FIS structure optimization method selection, and teaching signal selection. 

In one embodiment, the user makes the selection of fuzzy model, including one or 
more of: the number of input and/or output variables; the type of fuzzy inference model (e.g., 
Mamdani, Sugeno, Tsukamoto, etc.); and the preliminary type of membership functions. 

In one embodiment, a Genetic Algorithm (GA) is used to optimize linguistic variable 

25 parameters and the input-output training patterns. In one embodiment, a GA is used to 
optimize the rule base, using the fuzzy model, optimal linguistic variable parameters, and a 
teaching signal. 

One embodiment, includes optimization of the FIS structure by using a GA with a 
fitness function based on a response of a model of motorcycle and rider system. 
30 One embodiment, includes optimization of the FIS structure by a GA with a fitness 

function based on a response of the actual motorcycle/rider system. 
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The result is a specification of an FIS structure that specifies parameters of the optimal 
FC according to desired requirements. 

Brief Description of the Figures 
5 The above and other aspects, features, and advantages of the present invention will be 

more apparent from the following description thereof presented in connection with the 
following drawings. 

Figure 1 shows a robust intelligent control system of motorcycle and rider navigation 

model. 

10 Figure 2A is a photograph of an electric scooter used for demonstration purposes. 

Figure 2B shows a computer model of the scooter from Figure 2A 

Figure 3 shows a simulation model of the scooter from Figure 2B. 

Figure 4 shows a model of a spring and damper system. 

Figure 5 shows a model for aerodynamic forces (drag and lift). 
1 5 Figure 6 shows a tire model. 

Figure 7 shows a tire force coordinate system (Xt,Yt,Zt). 

Figures 8A-B show tire modeling using a JARI tire, including tire data (A), and tire 
model (B). 

Figure 9 shows an eigenvector of weave mode. 
20 Figure 10 shows a locus of eigenvalues of the weave and wobble mode for various 

velocities. 

Figure 1 1 shows features of the SC Optimizer. 
Figure 12A is a flowchart of the SC Optimizer. 
Figure 12B is a flowchart for using SC Optimizer. 
25 Figure 12C shows a flowchart for evaluating the SC Optimizer. 

Figure 13 shows a structure of the fuzzy control system and development of the fuzzy 
controller. 

Figure 14 shows a structure of the feedforward controller and its development. 
Figures 15A-B show geometrical parameters of a rider model and course, including 
30 reference roll angle (A), and length of reference and deviation (B). 
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Figure 16 shows the structure of a simulation to gain the teaching signals used by the 
SC optimizer. 

Figures 17A-G show examples of the data during fitness function evaluation 
simulation. 

5 Figures 18A-C show example results of the optimized parameters for different initial 

conditions, including Kpl (A), Kp2 (B), and Kd2 (C). 

Figure 19 shows membership functions for the fuzzy controller. 
Figure 20 shows outputs of the fuzzy controller. 
Figure 21 shows a simulation model for the fuzzy controller. 
10 Figures 22A-B show the map of the actual test course and simulated course. 

Figures 23A-F show the results of the simulation with the fuzzy controller (without 
excitation). 

Figures 24A-E show the results of the simulation with fuzzy controller (with 
excitation). 

15 Figures 25A-D show simulation results with different transport delay. 

Figure 26A-B show the look-forward model and the process to generate a feedforward 
rider model, where (A) shows the relation between deviation of course at length of reference 
and turning radius, and (B) shows the process of assembling the g feedforward rider model. 
Figure 27 shows a model for circular course simulation. 
20 Figure 28A-F show the optimized parameters. 

Figure 29A-B shows the relation between computed reference roll angle and the 
optimized parameter. 

Figure 30A-B show course lane change and model for lane change simulation. 
Figure 31 shows the feedforward model. 
25 Figure 32A-E show results of the simulation with a hold torque feedforward 

controller. 

Figure 33A-D show results of the simulation with transport delay. 
Figure 34 shows the structure of the measurement system. 
Figure 35 shows a test scooter with telemetry system. 
30 Figure 36A-B show an example of monitoring display (A) and measured data (B). 
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Detailed Description 

Figure 1 is a block diagram of an intelligent motorcycle (with rider) navigation control 
system based on soft computing and using a soft computing optimizer to generate a 
Knowledge Base (KB) 108 for a Fuzzy Controller (FC) 109. The motorcycle and rider system 
5 is represented by a block 102. Output from the block 102 is provided to a fitness function in 
block 104 of a GA in block 105. The objective information for the knowledge base (KB) 108 
is extracted (from the output of block 102) by a Genetic Analyzer (GA) 105 using the results 
of stochastic simulation of the motorcycle/rider dynamic behavior. The GA 105 and Fitness 
Function 104 are components of a System Simulation of Control Quality (SSCQ) 242. 

10 Using a set of inputs, and a fitness function in block 104, the GA in block 105 works in 

a manner similar to an evolutionary process to arrive at a solution, which is, hopefully, 
optimal. The GA in block 105 generates sets of "chromosomes" (that is, possible solutions) 
and then sorts the chromosomes by evaluating each solution using the fitness function in block 
104. The fitness function in block 104 determines where each solution ranks on a fitness scale. 

15 Chromosomes (solutions), which are more fit, are those, which correspond to solutions that 
rate high on the fitness scale. Chromosomes, which are less fit, are those, which correspond to 
solutions that rate low on the fitness scale. 

Chromosomes that are more fit are kept (survive) and chromosomes that are less fit are 
discarded (die). New chromosomes are created to replace the discarded chromosomes. The 

20 new chromosomes are created by crossing pieces of existing chromosomes and by introducing 
mutations. 

The motorcycle/rider model 102 is controlled by a linear P(D) controller 101. In one 
embodiment, the controller 100 uses both a PD controller 1012 and a P controller 1011 to 
control different aspects of the motorcycle system. The P(D) controller 101 has a linear 

25 transfer function and thus is based upon a linearized equation of motion for the control object 
in block 102. Prior art GAs used to program P(D) controllers typically use simple fitness 
function and thus do not solve the problem of poor controllability typically seen in 
linearization models. As is the case with most optimizers, the success or failure of the 
optimization often ultimately depends on the selection of the performance (fitness) function. 

30 Evaluating the motion characteristics of a nonlinear plant is often difficult, in part due to 

the lack of a general analysis method. Conventionally, when controlling a plant with nonlinear 
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motion characteristics, it is common to find certain equilibrium points of the plant and the 
motion characteristics of the plant are linearized in a vicinity near an equilibrium point. 
Control is then based on evaluating the pseudo (linearized) motion characteristics near the 
equilibrium point. 

5 Computation of optimal control based on soft computing includes the GA in block 105 

as the first step of global search for optimal solution on a fixed space of solutions. The GA 
searches for a set of control weights for the plant. Firstly the weight vector K = {k l ,... 9 k n }is 

used by a conventional P(D) controller in block 101 in the generation of a signal S (K) which 

is applied to the plant. The entropy S(5 (K)) associated with the behavior of the plant on this 

10 signal is assumed as a fitness function to be minimized. The GA is repeated several times with 
regular time intervals in order to produce a set of weight vectors. The vectors generated by the 
GA in block 105 are then provided to a FNN and the output of the FNN to a fuzzy controller 
in block 109. The output of the fuzzy controller in block 109 is a collection of gain schedules 
for the linear controllers in block 101 that control the motorcycle 102. 

15 So in summary, designing the intelligent control system includes two stages. Stage 1 

involves finding a teaching signal (also referred to as teaching patterns). Stage 2 involves 
approximation of optimal control. 

Finding teaching patterns (input-output pairs) of optimal control by using the GA are 
aspects of the SSCQ 242, based on the mathematical model of the controlled motorcycle and a 

20 physical criteria of minimum of entropy production rate. 

Approximation of the optimal control (from Stage 1) by the corresponding Fuzzy 
Controller (FC). This FC is referred to as the optimal Fuzzy Controller. 

The first stage is the acquisition of a robust teaching signal of optimal control without 
the loss of information. The output of first stage is the robust teaching signal, which contains 

25 information about the controlled object behavior and corresponding behavior of the control 
system. 

The second stage is the approximation of the teaching signal by building of a fuzzy 
inference system. The output of the second stage is the knowledge base (KB) 108 for the FC 
109. 
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The design of the optimal KB 108 includes specifying the optimal numbers of input- 
output membership functions, their optimal shapes and parameters, and a set of optimal fuzzy 
rules. 

In one embodiment for the Stage 2 realization, the optimal FC 109 is obtained using a 
5 FNN with the learning method based on the error back-propagation algorithm. The error back- 
propagation algorithm is based on the application of the gradient descent method to the 
structure of the FNN. The error is calculated as a difference between the desired output of the 
FNN and an actual output of the FNN. Then the error is "back propagated" through the layers 
of the FNN, and parameters of each neuron of each layer are modified towards the direction of 

10 the minimum of the propagated error. 

However, back-propagation algorithm has a few disadvantages. In order to apply back- 
propagation approach it is necessary to know the complete structure of the FNN prior to the 
training optimization. The back-propagation algorithm cannot be applied to a network with an 
unknown number of layers or an unknown number nodes. The back-propagation process 

15 cannot automatically modify the types of the membership functions. Usually, the initial state 
of the coefficients for back-propagation algorithm is set up randomly, and as a result, the 
back-propagation algorithm often finds only a "local" optimum close to the initial state. 

In spite of the above disadvantages, the error back-propagation algorithm is commonly 
used in practice. For example, the Adaptive Fuzzy Modeler (AFM) provided by 

20 STMicroelectronics Inc. creates Sugeno 0 order fuzzy inference systems using the error back- 
propagation algorithm. 

The algorithm used in the AFM is based on the following two steps. In the first step, a 
user specifies the parameters of a future FNN such as the numbers of inputs and outputs and 
the number of fuzzy sets for each input/output. The AFM "optimizes" the rule base, using a 

25 so-called "let the best rule win" (LBRW) technique. During this phase, the membership 
functions are fixed as uniformly distributed among the universe of discourse, and AFM 
calculates the firing strength of the each rule, eliminating the rules with zero firing strength, 
and adjusting centers of the consequents of the rules with nonzero firing strength. It is possible 
during optimization of the rule base specify the learning rate parameter, depending on the 

30 current problem. In AFM there is also an option to build the rule base manually. In this case, 
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user can specify the centroids of the input fuzzy sets, and then according to the specification, 
the system builds the rule base automatically. 

In the second step, the AFM builds the membership functions. The user can specify the 
shape factors of the input membership functions. Shape factors supported by AFM are: 
5 Gaussian, Isosceles Triangular, and Scalene Triangular. The user must also specify the type of 
fuzzy operation in Sugeno model: supported methods are Product and Minimum. 

After specifying the membership function shape and Sugeno inference method, the 
AFM starts optimization of the membership function shapes, using the structure of the rules, 
developed during Step 1. There are also some optional parameters to control optimization rate 
10 such as a target error and the number of iterations, the network should make. The termination 
condition on the optimization is reaching of the number of iterations, or when the error 
reaches its target value. 

AFM inherits the weakness of the back-propagation algorithm described above, and the 
same limitations. The user must specify types of membership functions, number of 
15 membership functions for each linguistic variable and so on. The rule number optimizer in the 
AFM is called before membership functions optimization, and as a result the system often 
becomes unstable during membership function optimization phase. 

However, some of the weakness of system such as the AFM can be ameliorated by 
optimizing the structure of the KB before training the KB. 
20 In Figure 1, the motorcycle (MC) and rider model (control algorithm) are shown as a 

simulation model 102. Figure 2 A shows an example of a real MC (in this cas, an electric 
scooter) and Figure 2B shows a model image of the scooter from Figure 2A. 

Figure 3 shows one embodimetn of a simulation motorcycle/rider model developed by 
using SimMechanics (a MATLAB add-on). Various rider models, including, for example, a 
25 fuzzy controller model and a hold torque feedforward model, can be used. 

In one embodiment, the models describe the stability of the MC during a straight line 
run. The model descriptions of the MC in this approache is provided by conventional straight- 
line modeling techniques. Such straight-line models typically do not include a steering model 
but do typically inlcude a passive rider model (e.g., a hands on and hands off model). 
30 In one embodiment, the model is similar to a model constructed using the propular 

DADS modeling program from LMS International. In this model, the forward velocity is 
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fixed. The rider steering models structure is similar to the fuzzy controller model, but the 
parameters of the rider model use fixed values which can be selected by trial and error. 

In one embodiment, the MC model is constructed using DADS, with a rider steering 
model based on a look-forward model. The parameters are fixed for each course. Lane change 
5 simulation and sine course simulation are shown. 

In one embodiment, the control behaviour of motorcycle and rider is described. The 
motorcycle model of the simulation is substantially simialr to the model described above. The 
rider model is a look-forward model. The parameters were fixec for each course. Lane change 
simulation is shown. 

10 In one embodiment, various driving maneuvres such as cornering and double lane 

change with a driver control algorithm are used with a motorcycle model using the ADAMS 
program available from MSC Software Corporation. The driver control model has two parts: a 
velocity control model and a steering control model. In this steering control model, the desired 
camber angle is estimated and depends on the driving maneuvre. The parameters for the 

1 5 controller can be fixed. 

In one embodiment, the stability control of a motorcycle/rider system based on fuzzy 
control in conjunction with a GA and an auto-tuning method is used. The model used in 
simulation need not be a true motorcycle model. For example, an inverted pendulum hinged to 
a rotating disk can be used as a representative model of the motorcycle 

20 One embodiment provides an autonomous two-wheeled vehicle and simulation results. 

The rider model can be divided into two parts: standing control and directional control. The 
control algorithm of standing stability is based on feedback control with gain parameters that 
depend on velocity. Directional control is the target roll angle, which can be determined by 
predicted course error. In one embodiment a spring-damper is connected between the steering 

25 bar and the actuator for steering to stabilize the motorcycle (which makes the control method 
relatively more complex). 

The effect on the perception of riding comfort issue for changing mass center position, 
total weight on the MC, the wheelbase, the front fork rake angle and the front wheel trail 
distance of the MC is discussed. The tire model is described in connectin with Figures 6-8. 



Figure 3 shows a simulation model of the scooter shown in Figures 2A-B. The model in 
the Figure 3 was made using SimMechanics which is an optional tool for kinematics analysis 
using the popular MATLAB computer program. 
Table 1 shows the model notations and parameter definitions. 
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Table 1 : Used notations 



Model notations 


r rf : Drive Torq(Nm) 


(x r ,y rt z r ) : directional vector of the tire rotational axis in gloval coordinate(m) 


r, : Steer Torq(Nm) 


(vjc ,vy ,vz ) : velocity vector of tire center point in gloval coordinate(m/s) 


r 4oW :steer torq to hold roll angle(Nm) 


(vjr,,vy ( ,vz,) : velocity vector of tire center point in tire force coordinate(m/s)(Fig.7) 


r, : lower body controlling Torq(Nm) 


(Fxa *Fra >Fzet)'f QTCe vector of tire at contact point in tire force coordinate(N)(Fig.7) 


r 2 : upper body controlling Torq(Nm) 


(T Xa ,0,7*2,, ):moment vector of tire at contact point in tire force coordinate(Nm) 


K pX ,K p2 ,K gains 


f r<j0 :static lateral force of tire at contact point in tire force coordinated) 


v : forward velocity m Is) 


7*^0 .-static overturning moment of tire at contact point in tire force coordinate(Nm) 


v ^ : reference velocity(m Is) 


T^oistatic aligning torque of tire at contact point in tire force coordinate(Nm) 


$ : roll angle(rarf ) 


T Xt :moment of tire excepting rolling moment 


^ , : reference roll angle(ra</ ) 

$rv/o • calculated reference roll angle(rad ) 

tb . ' correct inn fnr reference roll tmo\e( rnd \ 


around Xt(Fig.7) direction in tire force coordinate at center point of tire (Nm) 


T rt :moment of tire excepting rolling moment 


around Yt(Fig.7) direction in tire force coordinate at center point of tire (Nm) 


(F xt F r ,F Z ):force vector of tire at center point of tire in global coordinate(N) 


$ : rollrate(rarf Is ) 


(T x J r J 2 ):moment vector of tire excepting rolling moment 


a : required lateral acceleration(m/s 2 ) 


at center point of tire in global coordinate(Nm) 


g : gravity acceleration(m/s 2 ) 


M ^Tolling moment oftire(Nm) 


y : deviation at LR(m) 


K : Tire vertical spring constant(N/m) 


v 0 : deviation (course err)(m) 


D : Tire vertical damping(Ns/m) 


t : time(sec) 


extraction coeficient 


LR : length of reference(m) 


v, : forward speed for tirc(m/s) 


FF : fitness function 


a : coeficient for relaxation length(m) 


<o : yawrate(rarf Is) 


<r,:coef. for relaxation length proportion to load 


R : turning radius(m) 


ovcoef. for relaxation length constant 


9 : tire yaw angle in h@izontal plane(Fig.7) 


Cy ,,Cy M ,Cy l2 ,Cy 3 ,Cy 3M Cy 3J :coef. lateral force 


y : camber angle(rad)(Fig.7) 


Cz t ,Cz n ,Cz l2y Cz 3t Cz ilt Cz J2 :coef. overturning moment 


a : lateral slip(rad) 


Cx l9 Cx U9 Cx a , Cx 3t Cx u ,Cx„: coef. aligning torque 


k : rolling slip 


d F , t d F2l d Tlt d Tx : noises 


h : height of tire center(m)(Fig.7) 




tr : Tire original radius(m) 




te : Tire effective radius(m)(Fig.7) 




as : tire spin velocity(rad Is) 





The MC plant model includes: 

- 7 rigid bodies (mainframe in block 3001, rider in block 3002, front fork in block 
3003, front arm in block 3004, rear arm in block 3005, front & rear wheels in block 
3006 and in block 3007, respectively); 

- 5 rotating joints (rider roll in block 3008, steer in block 3009, rear arm pivot in 
block 3010, front & rear wheels in block 301 1 and in block 3012, respectively); 

- 1 slide joint (front suspension in block 3013); 

- 1 combination joint in block 3014 (a 6 degree of freedom main body to inertial 
system); 
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- 4 spring & damper (front suspension in block 3015, steer in block 3016, rider in 
block 3017, rear suspension in block 3018); 

- an external force in block 3019 (air force(drag & lift); and 

- a tire force in block 3020 for the front tire and a tire forces in block 3021 for the 
rear tire. 

Figure 4 shows the spring and damper model used in blocks 3015 and 3018. Parameters 
of the nonlinear spring, represented as look-up table, are provided to the block 4001. Damping 
coefficient are provided to the block 4002. 

Figure 5 shows the aerodynamic force model used in block 3019. Drag and lift 
parameters are provided to block 5001 and to block 5002, respectively. 

In the present disclosure, results from a MC simulation model are compared with results 
from an actual demonstration model. The differences between the MC model developed and 
the demonstrationmodel are : 

1. The linear tire force model in block 3020 and in block 3021 with first order lag 
was applied in the MC model (instead of the complex model used in the demonstration). 

2. Noise road signals were added to the tire model. 

3. Coordinate system was changed: +Z points up in the MC model and +Z points 
down in the demonstration. 

4. Front and rear brakes were removed in the MC model. 

5. Twist freedom of steering shaft was removed and only steering rotation 
freedom remains in block 3009. The removed twist axis was perpendicular axis to the 
steering shaft in central plane. The reason to remove the freedom is that in this model, 
linear tire model was used and small elastic deformation should be ignored in view of 
balance of accuracy. 

As mentioned above, blocks 3020 and 3021 represented the linear tire model. Tire 
properties are an important determinant of MC maneuverability and stability, and tire models 
are thus important to behavior simulations employing a mechanical analysis. Figure 6 is a 
block diagram of the tire model, where tire axis sensor data 6003 and tire spin sensor data are 
provided to a tire model 6001. The tire model 6001 provides rolling moment data 6005 and 
force moment and tire axis data 6006. The tire mathematical model in block 6001 is designed 
using the following equations: 
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9 = -arc tan(x r /y r ) (5-1) 
y = arcsin(z f /J(x r 2 +y 2 r +z 2 r ) ) (5-2) 
vx, = vx cos0+vy sinO (5-3) 
vy t = -vx sin#+vy cos# (5-4) 

vz, = vz (5-5) 

te =h / cosy (5-6) 

k - = -(v*, -te»m )lvx, (5-7) 

or = V jy, /vx, (5-8) 

=(K(tr-te)-Dte )(l+d Fz ) (5-9) 
^ =C,*F fc (5-10) 



o = (Cy x a +Cy iY ){\ +d Fy ) (5-11) 
T Zcl0 =(C Zl a+Cz 3 rXl+d T ,) (5-12) 
T Xa o=Cx l a+Cx 3 y(l+d Tx ) (5-13) 

(5-14) 
(5-15) 
(5-16) 
(5-17) 
(5-18) 
(5-19) 

(5-20) 

(5-21) 

(5-22) 

(5-23) 



Cy t =Cy u F Za0 +Cy l2 

c y i =c y i x F zc,o + cy n 

C Zl =Cz n F Za0 +Cz l2 

C Z 3~C Z 3 \F Zct 0~*~Cz n 

Cx y =Cx x x F Zct 0 +Cx 12 
Cx 3 =Cx 3l F Zct 0 +Cx 32 



— F + F = F 

T ' 1 Yet ~ 1 * ^ 



Yet 



YctO 



—T +7 = 7 

rr 1 ~* get 



y za 

—T + 7 = 7 

" 1 Xct ^ 1 1 



ZctO 



Xct 



XctO 
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M roil = T Zct smy-F Xct 4e (5-24) 

F x = F Xa cos 0-F Yct sin0 (5-25) 

F Y = F Xcs sin 6 + F Ya cos 0 (5-26) 

F z =F Zct (5-27) 
T x< = T x« + F Ya -tecosy + F Zct -tesiny (5-28) 

^ =-T Za cosy sin y (5-29) 

7\ = 7^, cos 0 -T Yt sin 0 (5-30) 

T Y =T Xt sir\0+T Yt cosd (5-31) 

T z = T Za cos 2 y (5-32) 



Inputs for the tire model in block 6001 are tire spin velocity in -block 6002 and the data 
of tire axis in block 6003. The data of the tire axis in block 6003 contain the height h of tire 

center and vector of tire axis (x r ,y r ,z r ) and velocity vector of tire center (v^v^,^). From 

5 these data, tire yaw angle, camber angle, velocity vector of tire center in tire force coordinate, 
tire effective radius, rolling slip, and lateral slip are calculated by Eqs (5-1) - (5-8) (see Figure 
7 for details). 

Static tire forces and moments at contact point in tire force coordinate are calculated by 

Eqs (5-9) - (5-20). Noise signals in block 6004 can be added to the tire model during these 
10 calculations. Output data from the tire model 6001 include: vertical force, driving force, static 

side force, static self aligning torque, static overturning torque at the contact point, and 

relaxation length. Conversion from static force and moments to dynamic moments is provided 

by Eqs (5-21) ~(5-23). 

The rolling moment in block 6005 for tire can be split by Eq. (5-24). The conversion 
15 from forces and moments at contact point to that at the shaft of tires in global coordinate 

system in block 6006 is provided by Eqs (5-25) - (5-32). 

The tire data for a desired tire can be measured or obtained, for example, from 

organizations such as the Japan Automobile Research Institute (JARI). In one embodiment, 

tires of the same specification can be used for the front and rear. 
20 In one embodiment, input data from a JARI tire model are vertical force, camber angle, 

and slip angle. Output data are side force, self aligning torque, overturning torque, and 

relaxation length. Figure 8A shows a sample tire data file containing data provided by JARI. 

This tire model is a linear model and thus less accurate than non-linear models for large 
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camber angles and large slip angles. Nevertheless, the linear model is adequate for many 
circumstances, including, but not limited to, modeling the scooter of Figure 2A in typical 
operating regimes. 

Figure 8B shows a transformation of the JARI tire model to add caculation of transient 
5 effects to the tire model. In Figure 8B input blocks, provide for arranging input data from 
sensed data of posistion, velocity, and posture of the tire. Additoinal blocks are provided to 
convert output data at tire contact in tire coordinate system to force and moment at tire center 
in inertia coordinate system containing transient effect. Forward and backward force 
calculation are provided to evaluate drag and driving forces. 

10 The model of the MC in SimMechanics is a complicated numerical simulation model 

that predicts the operation of the motorcycle system. In one embodiment, Eigenvalue analysis 
is used in connection with the numerical model to extract and explore various characteristics 
of the MC model produced by the simulation. 

For eigenvalue analysis, linearization of the model in the straight balanced running 

15 model can be used. The model has 12 degrees of freedom and adding differentiation to those 
12 degrees, gives are 24 degrees of freedom. Moreover, there are six degrees of freedom that 
originated in the feedback model of the relaxation length used in order to take the dynamic 
character of a tire into consideration. As a result of linearization of the model, a 30x30 
matrix representation of the system is obtained, correspdingly, 30 eigenvalues and 30 

20 eigenvectors are calculated. The characteristics of the motion of the model are determined by 
a relatively few modes among the 30 modes. The MC weave modes and wobble modes can be 
unstable, and the corredponding eigenvalues and eigenvectors are relatively large. 

Figure 9 shows the eigenvector of weave mode corresponding to 7 m/sec. The mode is a 
coupling mode that icludes cyclic motions of sway, roll, yaw and swing of steer. Yaw rate and 

25 swing of steer are in phase and roll angle is out of phase. 

Figure 10 shows the locus of eigenvalues of weave and wobble modes for various 
velocities. The weave mode is unstable for velocities under 10 m/sec. Wobble 1 mode 
corresponds to a simple swing of the steering. When the damping coefficient of steering is 
small, the mode will be unstable. In this model the damping coefficient C 5t = 5 is adopted. 

30 Wobble 2 mode includes rider roll, chassis roll, and steering swing. Wobble 3 mode includes 



yaw, roll, and steering swing. Table 2 shows eigenvalues of weave and wobble 1 when the 
damping coefficient of steering changes. 



Coefficient C tt 


Weave eigenvalue 


Wobble 1 eigenvalues 


0 


0.55 + 3.3i 


12 + 57i 


1 


0.63+3.3i 


8.0 + 58i 


3 


0.76 + 3.2i 


1.2 + 60i 


5 


0.88 + 3.0i 


-6.7 + 6H 



Table 2: 



5 

The results of simulation by the rider models were both stable. Robustness of the 
controllers was confirmed by simulation taking into acount the transmission delay of the 
response of actuator and noises in tire forces (as described below). 

In one embodiment, Stage 2 (defining the structure of the KB) is realized using a new 
10 structure based on Soft Computing (SC) techniques. The structure of the intelligent control 
system is based on the SC optimizer shown in Figure 1. In this case, the SC optimizer is used 
in place of the FNN block . 

The SC optimizer generates a KB of Fuzzy Inference System (FIS) from the digital in- 
out data. The SC optimizer provides GA based FNN learning, including rule extraction and 
15 KB optimization. In one embodiment, the SC optimizer can use as a teaching signal both 
pattern file and the model (control object) output. 

Figure 1 1 shows the features of the SC optimizer. The SC optimizer includes a fuzzy 
inference system (FIS) represented in the form of a FNN. In one embodiment, the SC 
optimizer allows structure selection, such as for example,Sugeno FIS order 0 and 1, Mamdani 
20 FIS, Tsukamoto FIS, etc. In one embodiment, the SC optimizer allows structure optimization 
selection, such as for example, Optimization via Genetic Algorithm (GA), Linguistic variables 
optimization, and rule-based optimization. In one embodiment, the SC optimizer allows 
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teaching source selection, such as for example: teaching signal as a look-up table (e.g., a table 
of in-out patterns); and/or teaching signal as a fitness function. When provided by a fitness 
function, the teaching signal fitenss function can be: calculated as a dynamic system response 
(e.g., using in MATLAB), calculated as a result of control of real control object (e.g., obtained 
5 using dSPACE or other method of connection of control object directly to computer), etc. 

Figure 12 is a flowchart for one embodiment of the SC optimizer. By way of 
explanation, and not by way of limitation, the operation of the flowchart divides operation in 
to five stages, shown as Stages 1, 2, 3, 4, and 5. 

In Stage 1, the user selects a fuzzy model by selection one of parameters such as, for 
10 example, the number of input and output variables, the type of fuzzy inference model 
(Mamdani, Sugeno, Tsukamoto, etc.), source of the teaching signal 

In Stage 2, a first GA (GA1) optimizes linguistic variable parameters, using the 
information obtained in Stage 1 about the general system configuration, and the input-output 
training patterns, obtained from the training signal as an input-output tables. 
15 In Stage 3 a precedent part of the rule base is created and rules are ranked according to 

their firing strength. Rules with high firing strength are kept, whereas weak rules with small 
firing strength are eliminated. 

In Stage 4, a second GA (GA2) optimizes a rule base, using the fuzzy model obtained 
in Stage 1, optimal linguistic variable parameters obtained in Stage 2, selected set of rules 
20 obtained in Stage 3 and the teaching signal. 

In Stage 5, the structure of FNN is further optimized. In order to reach the optimal 
structure, the classical derivative-based optimization procedures can be used, with a 
combination of initial conditions for back propagation, obtained from previous optimization 
stages. The result of Stage 5 is a specification of fuzzy inference structure that is optimal for 
25 the plant 120. Stage 5 is optional and can be bypassed. If Stage 5 is bypassed, then the FIS 
structure obtained with the GAs of Stages 2 and 4 is used. 

In one embodiment the teaching signal, representing one or more input signals and one 
or more output signals, is divided into input and output parts. Each of the parts is divided into 
one or more signals. Thus, in each time point of the teaching signal there is a correspondence 
30 between the input and output parts. 
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Each component of the teaching signal (input or output) is assigned to a corresponding 
linguistic variable, in order to explain the signal characteristics using linguistic terms. Each 
linguistic variable is described by some unknown number of membership functions, like 
"Large", "Medium", "Small", etc. 

"Vertical relations" represent the explicitness of the linguistic representation of the 
concrete signal, e.g. how the membership functions is related to the concrete linguistic 
variable. Increasing the number of vertical relations will increase the number of membership 
functions, and as a result will increase the correspondence between possible states of the 
original signal, and its linguistic representation. An infinite number of vertical relations would 
provide an exact correspondence between signal and its linguistic representation, because to 
each possible value of the signal would be assigned a membership function, but in this case 
the situations as "over learning" may occur. Smaller number of vertical relations will increase 
the robustness, since some small variations of the signal will not affect much the linguistic 
representation. The balance between robustness and precision is a very important moment in 
design of the intelligent systems, and usually this task is solved by Human expert. 

"Horizontal relations" represent the relationships between different linguistic variables. 
Selected horizontal relations can be used to form components of the linguistic rules. 

To define the "horizontal" and "vertical" relations mathematically, consider a teaching 
signal: 

[x(0,y(0], 

Where: 

t = 1,...,7V- time stamps; 

N - number of samples in the teaching signal; 

x{t) = ( jCj (f ), . . . , x m (t)) - input components; 

y (0 = 0>i (0> • • • y„ (0) - output components. 

Define the linguistic variables for each of the components. A linguistic variable is 
usually defined as a quintuple: (jc, T(x), U, G, M) , where x is the name of the variable, T(x) is 

a term set of the*, that is the set of the names of the linguistic values ofjc, with a fuzzy set 
defined in U as a value, G is a syntax rule for the generation of the names of the values of the 
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x and M is a semantic rule for the association of each value with its meaning. In the present 
case, x is associated with the signal name from xory , term set T(x) is defined using vertical 
relations, U is a signal range. In some cases one can use normalized teaching signals, then the 
range of U is [0,1] . The syntax rule G in the linguistic variable optimization can be omitted, 
5 and replaced by indexing of the corresponding variables and their fuzzy sets. 

Semantic rule M varies depending on the structure of the FIS, and on the choice of the 
fuzzy model. For the representation of all signals in the system, it is necessary to define m + n 
linguistic variables: 

Let[X,Y], X = (X l9 ...,X m ),Y = (Y l ,...J n ) be the set of the linguistic variables 

10 associated with the input and output signals correspondingly. Then for each linguistic variable 
one can define a certain number of fuzzy sets to represent the variable: 

Where 

fi J ^ , i = 1,..., in,./, = are membership functions of the i th component of the input 

15 variable; and 

fif f9 i = 1,..., nj i =l,...,/ y are membership functions of the z th component of the output 
variable. 

Usually, at this stage of the definition of the KB, the parameters of the fuzzy sets are 
unknown, and it may be difficult to judge how many membership functions are necessary to 
20 describe a signal. In this case, the number of membership functions (power of term set), 
l x € [1,L MAX ],/ = l,...,/w can be considered as one of the parameters for the GA (GA1) search, 

where L MAX is the maximum number of membership functions allowed. In one embodiment, 

Lmax is specified by the user prior to the optimization, based on considerations such as the 
computational capacity of the available hardware system. 
25 Knowing the number of membership functions, it is possible to introduce a constraint on 

the possibility of activation of each fuzzy set, denoted as p J X{ . 
One of the possible constraints can be introduced as: 

20 



p J x 'Z^,i = l 9 .. 9 m;j = l 9 ..J Xi 

This constraint will cluster the signal into the regions with equal probability, which is 
equal to division of the signal's histogram into curvilinear trapezoids of the same surface area. 
Supports of the fiizzy sets in this case are equal or greater to the base of the corresponding 
5 trapezoid. How much greater the support of the fuzzy set should be, can be defined from an 
overlap parameter. For example, the overlap parameter takes zero, when there is no overlap 
between two attached trapezoids. If it is greater than zero then there is some overlap. The 
areas with higher probability will have in this case "sharper" membership functions. Thus, the 
overlap parameter is another candidate for the GA1 search. The fuzzy sets obtained in this 

1 0 case will have uniform possibility of activation. 

Modal values of the fuzzy sets can be selected as points of the highest possibility, if the 
membership function has unsymmetrical shape, and as a middle of the corresponding 
trapezoid base in the case of symmetric shape. Thus one can set the type of the membership 
functions for each signal as a third parameter for the GA1 . 

15 The relation between the possibility of the fuzzy set and its membership function shape 

can also be found. The possibility of activation of each membership function is calculated as 
follows: 

A =p{^i=M j x i ) = j^f,M J xM(0) (6.1) 

Mutual possibility of activation of different membership functions can be defined from a 
20 geometrical view point as: 

where * denotes selected T-norm (Fuzzy AND) operation; 

j = l 9 ...J x J = 1,...,/^ are indexes of the corresponding membership functions. 

In fuzzy logic literature, T-norm, denoted as * is a two-place function from [0,l]x[0,l] 

25 to [0,1]. It represents a fuzzy intersection operation and can be interpreted as minimum 

operation, or algebraic product, or bounded product or drastic product. S-conorm, denoted 
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by+, is a two-place function, from [0,l]x[0,l] to[0,l]. It represents a fuzzy union operation 

and can be interpreted as algebraic sum, or bounded sum and drastic sum. Typical T-norm and 
S-conorm operators are presented in Table 3. 



T-norms (fuzzy intersection) 


S-conorms (frizzy union) 


min(x,y) - minimum operation 


max(;t, >/) - maximum operation 


xy- algebraic product 


x + y - xy - algebraic sum 


x * y = max [ 0, x + y- 1] - bounded 


;t + j; = min[l,;t + y]- bounded sum 


product 








x*y = < 


x, ify=\ 

y, if x=l - drastic product 
0, ifx,y <1 


x + y = < 


x, ify=0 

y, if x=0 - drastic sum 
0, if*,>;>0 



Table 3 



5 

If i = k , and j ± I , then equation (6.2) defines "vertical relations"; and if i * k , then 

equation (6.2) defines "horizontal relations". The measure of the "vertical" and of the 

"horizontal" relations is a mutual possibility of the occurrence of the membership functions, 

connected to the correspondent relation. 
10 The set of the linguistic variables is considered as optimal, when the total measure of 

"horizontal relations" is maximized, subject to the minimum of the "vertical relations". 

Hence, one can define a fitness function for the GA1 which will optimize the number 

and shape of membership functions as a maximum of the quantity, defined by equation (6.2), 

with minimum of the quantity, defined by equation (6.1). 
15 The chromosomes of the GA1 for optimization of linguistic variables according to 

Equations (6. 1) and (6.2) have the following structure: 

[^ i ,...j y J [^ i ,... > « r J[^ i ,...,r rw ] 

m+n m+n m+n 

Where: 

l X(Y)i e[l,L MAX ] are genes that code the number of membership functions for each 
20 linguistic variable X t {Y t ) ; 
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a X{Y)i are genes that code the overlap intervals between the membership functions of 
the corresponding linguistic variable X.ty) ; and 

T X(Y) are genes that code the types of the membership functions for the corresponding 
linguistic variables. 

Another approach to the fitness function calculation is based on the Shannon 
information entropy. In this case instead of the equations (6.1) and (6.2), for the fitness 
function representation one can use the following information quantity taken from the analogy 
with information theory: 

= ~P (*, I*, = M J x, ) l°g[p(xi |*, = M J 'x, )] ( 6 - 1 a) 



10 and 



= ~tM J x l (x i (t))\og[^ Xi (x i (t))] 

TV , =1 

x i\ x k \ ' x,=n> x .,x k =tt Xk J 

= -TfEIV (*,(')) * M' Xk W))] log[ M J Xl (x,(0) * M' Xt (** (0)1 



(6.2a) 



In this case, GA1 will maximize the quantity of mutual information (6.2a), subject to the 
minimum of the information about each signal (6.1a). In one embodiment the combination of 
information and probabilistic approach can also be used. 
15 In case of the optimization of number and shapes of membership functions in Sugeno - 

type FIS, it is enough to include into GA chromosomes only the input linguistic variables. The 
detailed fitness functions for the different types of fuzzy models will be presented in the 
following sections, since it is more related with the optimization of the structure of the rules. 
In one embodiment, a rules pre-selection algorithm selects the number of optimal rules 
20 and their premise structure prior optimization of the consequent part. 
Consider the structure of the first fuzzy rule of the rule base 
R l (t) = IF x x (t) is ri(x x ) AND x h (0 is M l 2 (x 2 ) AND-AND x m {t) is M } m (x m ) 
THEN y { (t) is /ifr^), y 2 (t) is fifcz\y 2 )> -.J'.CO « /& } 0O 
Where: 
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m is the number of inputs; 
n is the number of outputs; 
jc,(t),i = are input signals; 

yj(t) 9 j = are output signals; 

/4* are membership functions of linguistic variables; 
k = 1,.,.,/w + n are the indexes of linguistic variables; 

l k = 2,3,... are the numbers of the membership functions of each linguistic variable; 

//j'** - are membership functions of output linguistic variables, upper index; \l k } means 

the selection of one of the possible indexes; and 
t is a time stamp. 

Consider the antecedent part of the rule: 
R l m (t) = IF *i(0 is mI(x x ) AND x k (t) is ju\(x 2 ) AND- AND x m (t) is M l m (x m ) 
The firing strength of the rule R l in the moment t is calculated as follows: 

J4(0 = min[^ {x x (t)) 9 M l 2(x 2 (t)) 9 -,M l m (x m (t))] 
for the case of the min-max fuzzy inference, and as 

4(o=n[^^(o)^2(^(o),...,^(^(o)] 

for the case of product-max fuzzy inference. 

In general case, here can be used any of the T-norm operations. 

The total firing strength R l fs of the rule, the quantity R l fs (t) can be calculated as follows: 

for a continuous case, and: 

for a discrete case. 

In a similar manner the firing strength of each s -th rule is calculated as: 
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where 

m 

j = l,2,...,pj/ f is a linear rule index 

N - number of points in the teaching signal or maximum of / in continuous case. 
In one embodiment the local firing strength of the rule can be calculated in this case 
instead of integration, the maximum operation is taken in equation (6.3): 

R' fs =maxR' fs (t) (6.4) 

In this case, the total strength of all rules will be: 

R fs = Yj RS fs ' 

5=1 

where: 

m 

L 0 = Y[^k ' Number of rules in complete rule base 

Quantity R fs is important since it shows in a single value the integral characteristic of 

the rule base. This value can be used as a fitness function which optimizes the shape 
parameters of the membership functions of the input linguistic variables, and its maximum 
guaranties that antecedent part of the KB describes well the mutual behavior of the input 
signals. Note that this quantity coincides with the "horizontal relations," introduced in the 
previous section, thus it is optimized automatically by GAL 

Alternatively, if the structure of the input membership functions is already fixed, the 
quantities R* fi can be used for selection of the certain number of fuzzy rules. Many hardware 

implementations of FCs have limits that constrain, in one embodiment, the total possible 
number of rules. In this case, knowing the hardware limit L of a certain hardware 
implementation of the FC, the algorithm can select L < L 0 of rules according to a descending 

order of the quantities R* fs . Rules with zero firing strength can be omitted. 

It is generally advantageous to calculate the history of membership functions activation 
prior to the calculation of the rule firing strength, since the same fuzzy sets are participating in 
different rules. In order to reduce the total computational complexity, the membership 
function calculation is called in the moment t only if its argument jc(t) is within its support. 
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For Gaussian-type membership functions, support can be taken as the square root of the 
variance value <x 2 . 

In one embodiment, a rule pre-selection algorithm is based on a firing strength of the 
rules. The threshold level can be selected based on the maximum number of rules desired, 
based on user inputs, based on statistical data and/or based on other considerations. Rules with 
relatively high firing strength will be kept, and the remaining rules are eliminated. Rules with 
zero firing strength can be eliminated by default. In one embodiment, the presence of the rules 
with zero firing strength may indicate the explicitness of the linguistic variables (linguistic 
variables contain too many membership functions). The total number of the rules with zero 
firing strength can be reduced during membership functions construction of the input 
variables. This minimization is equal to the minimization of the "vertical relations." 

This algorithm produces an optimal configuration of the antecedent part of the rules 
prior to the optimization of the rules. Optimization of the consequential part of KB can be 
applied directly to the optimal rules only, without unnecessary calculations of the "un-optimal 
rules". This process can also be used to define a search space for the GA (GA2), which finds 
the output (consequential) part of the rule. 

A chromosome for the GA2 which specifies the structure of the output part of the rules 
can be defined as: 

[I, 1M ,=[/„...,/.],/, ={l,..J rk },k = l,...,n 

where: 

I, are groups of genes which code single rule; 

I k are indexes of the membership functions of the output variables; 

n is the number of outputs; and 
M is the number of rules. 

In one embodiment the history of the activation of the rules can be associated with the 
history of the activations of membership functions of output variables or with some intervals 
of the output signal in the Sugeno fuzzy inference case. Thus, it is possible to define which 
output membership functions can possibly be activated by the certain rule. This allows 
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reduction of the alphabet for the indexes of the output variable membership functions from 

to the exact definition of the search space of each rule: 

|imin /max \ . r imin /max rjmin /max ^ rimin /max ^ 

Thus the total search space of the GA is reduced. In cases where only one output 
5 membership function is activated by some rule, such a rule can be defined automatically, 
without GA2 optimization. 

In one embodiment in case of Sugeno 0 order FIS, instead of indexes of output 
membership functions, corresponding intervals of the output signals can be taken as a search 
space. 

10 For some combinations of the input-output pairs of the teaching signal, the same rules 

and the same membership functions are activated. Such combinations are uninteresting from 
the rule optimization view point, and hence can be removed from the teaching signal, reducing 
the number of input-output pairs, and as a result total number of calculations. The total 
number of points in the teaching signal (t) in this case will be equal to the number of rules plus 

15 the number of conflicting points (points when the same inputs result in different output 
values). 

The previous section described optimization of the FIS, without the details into the type 
of FIS selection. In one embodiment, the fitness function used in the GA2 depends, at least in 
part, on the type of the optimized FIS. Examples of fitness functions for the Mamdani, Sugeno 
20 and/or Tsukamoto FIS models are described herein. One of ordinary skill in the art will 
recognize that other fuzzy models can be used as well. 

Define error E p as a difference between the output part of teaching signal and the FIS 
output as: 

E^Ud'-FW^,.^)? and E = ^E" , 

p 

25 where xf 9 x p 9 ... 9 xfl and d p are values of input and output variables in the p training pair, 
respectively. The function F(xf 9 x p 9 ... 9 x p ) is defined according to the chosen FIS model. 
For the Mamdani model, the function F(xf , x{ , x p n ) is defined as: 
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zim<*> z** 

/=1 i=l /=l 

n 

where z' = !"[/*;,(*,■) m ^ y l is the point of maximum value (called also as a central value) 

1-1 

°f mJ 00 > 0 denotes * e selected T-norm operation. 

Typical rules in the Sugeno fuzzy model can be expressed as follows: 
IF x x is fPjfr) AND x 2 is fi (l) h (x 2 ) AND ... AND jc/is (jcJ 

TEENy = f l (x X9 ...,x n ) 9 

where / = 1,2,..., M - the number of fuzzy rules M defined as{ number of membership 
functions of x { input variable} x{ number of membership functions of x 2 input 

variable} x . . . x { number of membership functions of x n input variable} . 
The output of Sugeno FIS is calculated as follows: 

Fl*,,*, O'^ir^r 1 • (6-6) 

ZIKw 

l=\ i=\ 

Typical rules in the first-order Sugeno fuzzy model can be expressed as follows: 
IF is M il) jS x i) AND x 2 is M (l) j 2 ( x i) AND AND *■ is t* i0 Ja M 
THEN;; = f l (x X9 ...,x tt ) = p«\ + p { 2 l) x 2 + ...rf>x. + , 

(Output variables described by some polynomial functions.) 
The output of Sugeno FIS is calculated according equation (6.6). 

Typical rules in the zero-order Sugeno FIS can be expressed as follows: 
IF x x is y^Oc,) AND x 2 is fi {l) h (x 2 ) AND ... AND x n is H {1) Jn (x n ) 

THENj/ = r (/) , 

The output of zero-order Sugeno FIS is calculated as follows 
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l=\ 1=1 



Fl Xl , Xl ,..., Xm ).M-M (6.7) 



/ = 1 1 = 1 

The typical rule in the Tsukamoto FIS is: 
IF jc, is fi {i \ (x } ) AND x 2 is n {l) h (x 2 ) AND ... AND x n is /Jp\{x m ) 

THEN y is Ml l) (y)> 

5 where j\ e I is the set of membership functions describing linguistic values of x x 

input variable; j 2 e I is the set of membership functions describing linguistic values of x 2 
input variable; and so on , j n e I is the set of membership functions describing linguistic 
values of x n input variable; and keO is the set of monotonic membership functions 
describing linguistic values of y output variable. 
10 The output of the Tsukamoto FIS is calculated as follows: 



M n M 

Fix x \ - 1=1 1=1 - /=1 



Ul/ii*,) Y z ' (6-8) 

l=\ i=\ l=\ 

where z'=fl//> f ) ™ d z ' =^{y') 

1=1 

Stage 4 described above generates a KB with required robustness and performance for 
15 many practical control system design applications. If performance of the KB generated in 
Stage 4 is, for some reasons, insufficient, then the KB refinement algorithm of Stage 5 can be 
applied. 

In one embodiment, the Stage 5 refinement process of the KB structure is realized as 
another GA (GA3), with the search space from the parameters of the linguistic variables. In 
20 one embodiment the chromosome of GA3 can have the following structure: 

|[A 19 A 2 , A 3 ]j^ ; A f - e[-prmj ,1- prmj];i = l,2,3;y = 1,2,...,Z , where L is the total number of 

the membership functions in the system 
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In this case the quantities A, are modifiers of the parameters of the corresponding fuzzy 

set, and the GA3 finds these modifiers according to the fitness function as a minimum of the 
fuzzy inference error. In such an embodiment, the refined KB has the parameters of the 
membership functions obtained from the original KB parameters by adding the modifiers 
5 prm new i = prm.+A,.. 

Different fuzzy membership function can have the same number of parameters, for 
example Gaussian membership functions have two parameters, as a modal value and variance. 
Iso-scalene triangular membership functions also have two parameters. In this case, it is 
advantageous to introduce classification of the membership functions regarding the number of 

10 parameters, and to introduce to GA3 the possibility to modify not only parameters of the 
membership functions, but also the type of the membership functions, form the same class. 

GA3 improves fuzzy inference quality in terms of the approximation error, but may 
cause over learning, making the KB too sensitive to the input. In one embodiment a fitness 
function for rule base optimization is used. In one embodiment, an information-based fitness 

15 function is used. In another embodiment the fitness function used for membership function 
optimization in GA1 is used. To reduce the search space, the refinement algorithm can be 
applied only to some selected parameters of the KB. In one embodiment refinement algorithm 
can be applied to selected linguistic variables only. 

The structure realizing evaluation procedure of GA2 or GA3 is shown in Figure 12B. In 

20 Figure 12B, the SC optimizer 12001 sends the KB structure presented in the current 
chromosome of GA2 or of GA3 to FC 12101. An input part of the teaching signal 12102 is 
provided to the input of the FC 12101. The output part of the teaching signal is provided to the 
positive input of adder 12103. An output of the FC 12101 is provided to the negative input of 
adder 12103. The output of adder 12103 is provided to the evaluation function calculation 

25 block 12104. Output of evaluation function calculation block 12104 is provided to a fitness 
function input of the SC optimizer 12001, where an evaluation value is assigned to the current 
chromosome. 

In one embodiment evaluation function calculation block 12104 calculates 
approximation error as a weighted sum of the outputs of the adder 12103. 
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In one embodiment evaluation function calculation block 12104 calculates the 
information entropy of the normalized approximation error. 

In one embodiment of Stages 4 and 5 the fitness function of GA can be represented as 
some external function Fitness = f (KB) , which accepts as a parameter the KB and as output 
5 provides KB performance. In one embodiment, the function / includes the model of an actual 
plant controlled by the system with FC. In this embodiment, the plant model in addition to 
plant dynamics provides for the evaluation function. 

In one embodiment function / might be an actual plant controlled by an adaptive 
LINEAR controller with coefficient gains scheduled by FC and measurement system provides 
10 as an output some performance index of the KB. 

In one embodiment the output of the plant provides data for calculation of the entropy 
production rate of the plant and of the control system while the plant is controlled by the FC 
with the structure from the KB. 

In one embodiment, the evaluation function is not necessarily related to the mechanical 
15 characteristics of the motion of the plant (such as, for example, in one embodiment control 
error) but it may reflect requirements from the other viewpoints such as, for example, entropy 
produced by the system, or harshness and or bad feelings of the operator expressed in terms of 
the frequency characteristics of the plant dynamic motion and so on. 

Figure 12C shows one embodiment the structure-realizing KB evaluation system based 
20 on plant dynamics. In Figure 12C, and SC optimizer 12201 provides the KB structure 
presented in the current chromosome of the GA2 or of the GA3 to an FC 12301. the FC is 
embedded into the KB evaluation system based on plant dynamics 12300. The KB evaluation 
system based on plant dynamics 12300 includes the FC 12301, an adaptive LINEAR 
controller 12302 which uses the FC 12301 as a scheduler of the coefficient gains, a plant 
25 12303, a stochastic excitation generation system 12304, a measurement system 12305, an 
adder 12306, and an evaluation function calculation block 12307. An output of the linear 
controller 12302 is provided as a control force to the plant 12303 and as a first input to the 
evaluation function calculation block 12307. Output of the excitation generation system 12304 
is provided to the Plant 12303 to simulate an operational environment. An output of the Plant 
30 12303 is provided to the measurement system 12305. An output of the measurement system 
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12305 is provided to the negative input of the adder 12306 and together with the reference 
input Xref forms in adder 12306 control error which is provided as an input to the linear 
controller 12302 and to the FC 12301. An output of the measurement system 12305 is 
provided as a second input of the evaluation function calculation block 12307. The evaluation 
5 function calculation block 12307 forms the evaluation function of the KB and provides it to 
the fitness function input of SC optimizer 12201. Fitness function block of SC optimizer 
12201 ranks the evaluation value of the KB presented in the current chromosome into the 
fitness scale according to the current parameters of the GA2 or of the GA3. 

In one embodiment, the evaluation function calculation block 12307 forms evaluation 
10 function as a minimum of the entropy production rate of the plant 12303 and of the linear 
controller 12302. 

In one embodiment, the evaluation function calculation block 12307 applies Fast Fourier 
Transformation on one or more outputs of the measurement system 12305, to extract one or 
more frequency characteristics of the plant output for the evaluation. 

15 In one embodiment, the KB evaluation system based on plant dynamics 12300 uses a 

nonlinear model of the plant 12303. In one embodiment, the KB evaluation system based on 
plant dynamics 12300 is realized as an actual plant with one or more parameters controlled by 
the adaptive linear controller 12302 with control gains scheduled by the FC 12301. In one 
embodiment plant 12303 is a stable plant. In one embodiment plant 12303 is an unstable plant. 

20 The output of the SC optimizer 12201 is an optimal KB 12202. 

In Figure 13, the structure of the fuzzy controller system and its development stages is 
described, combining the following Stages: Stage 1 (13100), teaching signal acquisition stage 
for several initial conditions; Stage 2 (13200), arbitrary course simulation. 

In stage 1(13100), several combinations of reference velocities and roll angles are set, 

25 which are also used as initial velocities and roll angles. The plant model of motorcycle 13101 
(as shown in Figure 3) and PD controller for steering 13103 and P controller for driving 13104 
are embedded into the structure of Stage 1 (13100). The calculation for controllers 13103, 
13104 is based on the equation (7-l)-(7-2). And the gain parameters for the PD controller for 
steering 13103 and the P controller for driving 13104 are provided by the GA 13102. The GA 

30 13102 provides optimization of the parameters executed through the evaluation of the value of 
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the Fitness Function 13105 for various conditions. The calculation of the Fitness Function 
13105 is based on the equation (7-6). 

Optimized parameters are stored in the KB 13002. Data for the KB 13002 are 
transformed into fuzzy rule form for the KB 13003 by the SCO 13001. The same plant model 
5 of MC 13201 and PD controller for steering 13203 and P controller for driving 13204 are 
embedded into the structure of Stage 2 (13200). The calculation for the controllers 13203, 
13204 is based on equations (7-l)-(7-2). In this stage, the gain parameters for the PD 
controller for steering 13203 and P controller for driving 13204 are provided by the FC 13206 
using the KB of the fuzzy rule 13003. Input data for the P controller for driving 13204 is the 
10 difference between forward velocity and reference velocity. Input data for the PD controller 
for steering 13203 are the difference between roll angle and reference roll angle, which is 
calculated by cp re f calculator 13207. Input data for the calculator 13207 are LR 13208, course 
data 13212 and position and yaw of the plant model of the MC 13201. The calculation is 
based on equations (7-3)-(7-5). 

^=-K pl (v-v ref ) (7-1) 

15 r,=-tf, 2 (*-^)-*, 2 rf (7-2) 

t=LR/v (7-3) 
a=2y/t 2 (7-4) 

(f> ref =t an- 1 (a/g) = tan- 1 (^^) (7-5) 

2 

FF = (v -v rg/ ) 2 x0.2+(^-^ e/ ) 2 +^ +r/x0.0001+r s 2 x0.0005 (7-6) 

Figure 14 is a block diagram showing the structure of the hold steering torque 
feedforward system and its development. Figure 14 shows the following 3 stages. Stage 1 
20 (14100) is a data acquisition stage for several conditions. Stage 2 (14200) is a parameter 
optimization stage for several forward velocities. Stage 3 (14300) is an arbitrary course 
simulation. 

In Stage 1 (14100), several reference velocities are set. The reference velocities are also 
used as initial velocities. Several radius of circle course 14112 and several LR 14108 are used. 
25 As combination of several reference velocities and several radius of circle course 14112 and 
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several LR 14108, many cases of simulation must be executed. The plant model of the 
motorcycle 14101, the PD controller for steering 14103, and the P controller for driving 14104 
are embedded into the structure of Stage 1 (14100). The calculation for controllers 14103, 
14104 is based on equation (7-l)-(7-2). The P gain parameter for PD controller for steering 
5 14103 and P controller for driving 14104 are fixed. The D parameter for PD controller for 
steering 14103 are provided by a GA 14102. The GA 14102 feeds the value of <p co i 14109 and 
ihoid 14110. (pcoi 14109 is used as input for q> re f ca ic- 14107. Other input data of the calculator 
14107 are LR 14108 and course data (a circular course is used here) 14112 and position and 
yaw of the plant model of the motorcycle 14101. The calculation 14107 is based on equations 

10 (7-7)-(7-lll). The value Thoid 14110 is added to the output of the PD controller for steering 
14103. By using the GA 14102, the optimization of the parameters is executed through the 
evaluation of the value of the Fitness Function 14105 for each desired condition. The 
calculation of Fitness Function 14105 is based on equation (7-12). Optimized parameters of 
cpcoi 14109 and x ho id 141 10 are stored in the KB 14001. 

15 The data of KB 14001 are used in the Look up table 14211 in Stage 2 (14200). The 

Look up table 14211 brings optimized (p co i 14209 and T ho id 14210 according to LR 14208 and 
reference roll angle calculated by <p re f caic- 14107 and forward velocity of the plant model of the 
motorcycle 14201. In stage 2 (14200), several reference velocities are set, which are also used 
as initial velocities. The plant model of the motorcycle 14201, the PD controller for steering 

20 14203, and the P controller for driving 14204 are embedded into the structure of Stage 2 
(14200). The calculation for controllers 14203, 14204 is based on equations (7-l)-(7-2). The 
gain parameters for P controller for driving 14204 are fixed. In this stage, the gain parameters 
for the PD controller for steering 14203 are provided by a GA 14202. The LR 14208 is also 
provided by the GA 14202. The LR 14208, course data (lane change course which is shown in 

25 Fig.30(a) is used here) 14212 and position and yaw of the plant model of motorcycle 14201 
and optimized (p co i 14209 are used as inputs for (p re f ca ic 14207. The calculation 14207 is based 
on equations (7-7)-(7-l 1). The value ihoid 14210 is added to the output of the PD controller for 
steering 14203. By using the GA 14202, the optimization of the parameters is executed 
through the evaluation of the value of the Fitness Function 14205 for each reference velocity. 

30 The calculation of the Fitness Function 14205 is based on equations (7-13). Optimized 
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parameters of the gain parameters for the PD controller for steering 14203 and LR 14208 are 
added to the KB 14001, and the KB 14002 is constructed. 

The data of the KB 14002 are used in the Look-up table 14311 in Stage 3 (14300). The 
Look-up table 14311 brings optimized q> co i 14309 and Thoid 14310 and the gain parameters for 
5 the PD controller for steering 14303 and LR 14308 according to reference roll angle 
calculated by q> re fcaic 14107 and forward velocity of the plant model of motorcycle 14301. The 
plant model of the motorcycle 14301, the PD controller for steering 14303, and the P 
controller for driving 14304 are embedded into the structure of Stage 3 (14300). The 
calculation for controllers 14303, 14304 is based on equations (7-1) -(7-2). 

10 The gain parameters for the P controller for driving 14304 are fixed. In this stage, the 

gain parameters for the PD controller for steering 14303 are provided by the Look-up table 
14311 as mentioned above. The LR 14308, arbitrary course data 14312, and the position and 
yaw of the plant model of the motorcycle 14301 and optimized (p co i 14309 are used as inputs 
for (p re f caic 14307. The calculation 14307 is based on equation (7-7)-(7-ll). The value x ho id 

15 14310 is added to the output of the PD controller for steering 14303. 



o)=v/R (7-7) 
a=v*6) (7-8) 

R =(y 2 +LR 2 )/2y (7-9) 

<p refQ =tan" 1 (a/g) = tan- 1 ( H (7-10) 

wo sj g *(y 2 +LR ) 

kef =<f>ref0-<l>cor O A I) 

FF = y 0 2 +(<f>-<f> ref ) 2 (7-12) 
FF =y 0 2 +r 5 2 x0.007 (7-13) 



20 The optimized KB defines the optimal fuzzy control sytem to realize a robust rider 

model, which can drive the motorcycle (plant model) passing through an arbitrary course trace 
with arbitrary velocity without fall down, and within the desired control error tolerance and 
control quality. 
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In one embodiment, the rider model presented is a combination model of look-forward 
a model and an ideal roll control model. To calculate steering torque, using deviation at 
length of reference is from look-forward model and using roll rate is from ideal roll control 
model (see Figure 15). The rider model uses a fuzzy controller that optimizes the gain 
parameters corresponding to the velocity and the roll angle. To optimize the gain parameters, 
the Soft Computing Optimizer based on GA (genetic algorithm) with different types of fitness 
function is used in a simulation to control the velocity and roll angle to reference value. The 
method using the Soft Computing Optimizer is divided into two parts. The first part is to gain 
the teaching signals by simulation, and the second part is to make the optimized fuzzy 
controller by using the teaching signals. 

The equations for controller torque are: 

*,«-*>-v^) (9-1) 

(9-2) 

The equations for reference roll angle are: 

t=LRIv (9-3) 
a=2y It 1 (9-4) 

<fi ref =tan- 1 (a/^) = tan- I (^^l) (9-5) 

The equations for steer torque and body movement controller are: 

r s =K^ + K ys y (9-6) 
t^Kj, (9-7) 
r 2 =K^ + K y2 y (9-8) 

The reference roll angle is given by course radius: 

r d -K pV (v ref -v) (9-9) 

To obtain the teaching signal, the Soft Computing Optimizer uses a program mexGA, 
a MATLAB program script.m a MATLAB program fitness.m and the MC model. The 
structure is shown in Figure 16. In script.m the initial conditions and reference values are set. 
In fitness.m the fitness function is defined as: 

M-file: file in MATLAB language 

Reference velocity = 5 7 10 15 m/s 

Reference roll angle= 0° 10° 30° 
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Optimized gain Kp 1 ,Kp2,Kd2 
Fitness function 

2 

FF =(v -v ref ) 2 x0.2+(<f>-<f> ref ) 2 +<f> +r/ x 0.000 1+r/x 0.0005 (9-11) 

GA parameter 
Population = 100 
Generation = 5 
Mutation_rate = 0.6 
Crossover_rate = 0.9 
Delete_rate =0.8 
Fitness_reduce =0.05 

Figure 17 shows examples of teaching signals obtained from the simulation. Figure 18 
shows results of optimized parameters for different initial conditions. Several features of the 
teaching signals produced by simulation are noted as follows. Kpl is relatively large when the 
15 velocity is fast, regardless of roll angle. Kp2 change relatively little except when velocity is 
slow and roll angle is large. When the reference roll angle is 30 deg. and the reference velocity 
is 5m/s, the resulting roll angle converges at less than 27deg., regardless of gain parameters. 
This appears to be the basis for the preceding exception. Torque of the front tire around the 
steering axis is largely the reason for this phenomenon. Kd2 is relatively large when the 
20 velocity is slow, regardless of roll angle. 

As a next step, the SC optimizer computes the optimum structure of the KB for the 
fuzzy controller. The input for the SC optimizer is the teaching signal. The membership 
functions for the KB are set to make the input data of fuzzy operation from the input data of 
the fuzzy controller, such as velocity and roll angle. The SC optimizer determines functions 
25 and makes fuzzy rules and optimizes them by the teaching signals using GA. Figures 19 and 
20 shows the results of membership function design and FC output, correspondingly. 

In Figure 21 is a block diagram of the simulation model including a fuzzy controller 
block 2101 The block 2102 serves as the external performance function of the SC optimizer, 
which provides gain parameters for the simulation using the fuzzy rules. A reference roll angle 
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calculator 2103 block calculates the reference roll angle from assigned course data and LR 
(length of reference) and input data(position, velocity and yaw). 

The test course assigned here and the simulation are shown in Figure 22. The course is 
driven in a counterclockwise direction. The course starts in a straight line and goes into a J- 
5 turn of 13 meter radius. After a short straight and turn, it goes to lane change and J-turn of 10 
meter radius. The course then goes through a large radius turn and back the starting straight. 
In this simulation, the motorcycle run along this course two times, and as the reference 
velocity the experiment data (Figure 23A) was used. In the first turn the reference velocity is 
approximately 5m/s and in the second turn the reference velocity is approximately 7m/s. LR is 
10 6 meters for the first turn and is 7.5 meters for the second turn. 

Figures 23B-E show the resulting data of simulation and corresponding measured 
data. At the 10 meter J-turn RIO in the second circle (in time scale 70sec), the simulation data 
vibrates strongly compared to experiment data. The reason lies in the torque of the front tire 
around the steering axis. When it is about to go into a regular circle, the calculated control 
15 torque approaches 0. If the torque of the front tire around the steering axis is 0, the regular 
circle revolution would be kept. But the torque is not 0, so the vibration is provoked. Fig. 23F 
shows the change of gain parameters by the fuzzy controller. 

To evaluate the robustness of the controller, a simulation with disturbance in tire load 
was carried out as shown in Figures 24A-E. The effect is relatively small as the frequency of 
20 the weave mode, which is the dominant mode of MC at researched velocity range, is about 
0.5Hz. So the high frequency noise effect is small. If the velocity range was different and the 
high frequency modes were dominant, the influence of the tire noise would be more 
important. 

The other tests for the robustness was simulations assuming the transport delay of the 
25 actuator. As shown in Figures 25A-D, until 0.04 sec of delay is introduced, the result did not 
change much. However, when assuming a 0.05 sec delay, the vibration of the MC diverged 
and the simulation became unstable. 

In optimization of control of a motorcycle, stabilizing roll angle goes together with 
steering the desired course (see Figure 26A). It is thus useful to first explore control methods 
30 that can stabilize the motorcycle when driven in a circle. The reason is that roll angle will be 
almost constant in a circular path for constant velocity and any course desired can be resolved 



into combinations of short circular arcs. If the control method can stabilize the motorcycle 
through a circular arc of any desired radius, then the MC can be stably steered through any 
course constructed from arcs having radii within the design range. 

Accordingly, a circular course was used for the simulation and the circular course was 
run while computing a reference roll angle, and optimization of the parameter by the SC 
optimizer was used during the simulation (see Figure 26B). As described by equation 9-4, it is 
necessary to provide a holding torque to counter the torque of the front tire around the steering 
axis to hold the reference roll angle. So, a feedforward rider model can be used. The method 
of computing the reference roll angle in equations (9-3)-(9-5) assumes the course deviation to 
be small. Ignoring the change of yaw, the constant lateral acceleration cancels the predicted 
course deviation. So when the course deviation is large, computing reference roll angle by this 
method has a large error. The present method for computing reference roll angle (equations 
(9-14)-(9-17) assumes the arc course. So this method fits the circle course rather better. 
However, there is still some difference between computed reference roll angle and actual 
converged roll angle. The correction value introduced in equation (9-18) and optimized in 
simulation of circle course reduces the error. The correction compensates for the overturning 
moment of the tire and wheel base length. The optimization of parameter (the holding torque 
and the correction value of reference roll angle) by simulation of circle course, the 
controllability for transient response can not be optimized enough. The lane change simulation 
was added for optimization of length of reference(LR) and the PD gains(Kp2,Kd2). 

In this method velocity control gain was fixed. The equations for controller torque are: 



The equations for reference roll angle are: 



(9-13) 



(9-12) 



co =v IR 
a -v * co 
from fi g. 26(a) 
R = (y 2 +LR 2 )/2y 



(9-14) 
(9-15) 



(9-16) 



^ /0 =tan- , (a/^) = tan- , C 



2y *v 2 



(9-17) 



8 



*(y 2 +LR 2 ) 



<l>ref =<f>refO-</>c 



(9-18) 



39 



Figure 27 shows the model for circular course simulation. 
To construct the feedforward rider model, simulation of the circular course controlled 
by SC optimizer was done under several conditions. The holding torque and the correction 
value of the reference roll angle were optimized and the converged roll angle were obtained in 
5 each condition. Gain parameters Kpl and Kp2 were fixed (Kp 1=1 000, Kp2=80). The gain 
parameter Kd2 was added to the optimization. 

The results of the simulation are shown in Figure 28. The holding torque and the 
correction value of the reference roll angle and the converged roll angle become larger as the 
radius of the course becomes smaller. LR affects the correction value of reference roll angle. 
10 The holding torque and the correction value of reference roll angle are nearly in proportion to 
the computed reference roll angle (Figure 29A-B). 

It is difficult to optimize the controllability of transient response by circle course. The 
lane change simulation (see Figure 30) can be used for parameters to be optimized in transient 
conditions. The result of the optimization in circle course simulation was used in lane change 
15 simulation. 

ref. velocity , / 

optimized parameter Kp2 Kd2 

fitness function 

20 FF =>> 0 2 +r f 2 x0.007 (9-19) 

Holding torque was computed by interpolation of input data as a function of the 
computed reference roll angle, (assuming velocity as constant) Correction of reference roll 
angle was computed by interpolation of input matrix data as function of LR and computed 
reference roll angle, (assuming velocity as constant) 
25 The feedforward rider model (see Figure 31) was made as the table data (Table2) 

which are optimized by simulations of the circle course and the lane change course controlled 
by the SC optimizer. The resulting data of simulation in the test course fit well to experiment 
data (Figures 32A-E). At the 10 meter J-turn in the second circuit (in time scale 70sec), the 
simulation data does not vibrate. 



The tests for the robustness in the simulation are provided by assuming the transport 
delay of the steering actuator. As shown in Figures 33A-D, until 0.05 sec of delay is 
introduced, the control result does not change substantially. Introducing a 0.07 sec delay 
causes vibration of the MC and instability in the simulation 

Figure 34 shows the sensor and measurement apparatus. The measuring system has 
two portions. The first portion 34200 is integrated into the test scooter as shown in Figure 35. 
The second portion 34100 is the measurement base which includes a wireless modem 34101, 
GPS reference station 34102, transmitter board 34103, and computer with monitor 34104. The 
measured data (position and attitude of the scooter, steering angle, steering torque) can be 
seen on the monitor as shown in Figure 36. The first portion 34200 integrated in the test 
scooter includes a wireless modem 34201, steering sensor system 34202, which is installed in 
the steering part of the scooter as shown in Figure 35, driving motor sensor 34203, position 
and attitude sensors 34204 attached to the scooter as shown in Figure 35, and ECU unit 
34205. The block 34202 and 34203 are the additional part to the original system in order to 
specialize the original system for the scooter measurement. 
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